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ABSTRACT 


This report describes the planning, observations and 
adjustment of high-precision gravity survey networks established 
on the islands of Maui and Oahu as part of the geodetic- 
geophysical program in support of lunar laser ranging at 
Haleakala, Maui, Hawaii. 

The gravity survey networks include 43 independently 
measured gravity differences along the gravity calibration 
line from Kahului Airport to the summit of Mt. Haleakala, 
together with some key points close to tidal gauges on Maui, 
and 40 gravity differences within metropolitan Honolulu on 
Oahu . 


The results of our 1976-1978 survey are compared with 
surveys made in 1961 by the 1381st Geodetic Survey Squadron, 

Air Photographic and Charting Service, United States Air 
Force and with those made in 1964-1965 by the personnel of 
Hawaii Institute of Geophysics on the islands of Maui and Oahu. 

All final gravity values are given in the system of the 
International Gravity Standardization Net 1971 (IGSN 71); 
values are obtained by subtracting 14.57 mgal from the Potsdam 
value at the gravity base station at the Hickam Air Force Base, 
Honolulu. For the annual comparison of inter-island gravity 
surveys the secondary gravity base station at the Kahului 
Airport, Maui has been selected as the reference station. The 
tie between the two gravity networks was made from Kahului 
Airport to the inter-island terminal at Honolulu International 
Airport . 
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I. INTRODUCTION 


To verify and measure any tectonic plate motion between 
the Hawaiian Islands and the western United States mainland 
and to collect data for polar motion studies, the United States 
lunar laser ranging program maintains two permanent lunar laser 
ranging installations. The McDonald Observatory, on the North 
American tectonic plate, is located at the summit of Mt . Locke 
in the Davis Mountains of western Texas. It is utilized by 
the Marine Science Institute, University of Texas. LURE 
Observatory, on the Pacific plate, is located at the summit 
of Mt . Haleakala on the Island of Maui. It is utilized by 
the Institute for Astronomy, University of Hawaii. 

To obtain maximum information about plate motion from 
lunar laser ranging in a volcanically active region such as 
the Hawaiian Islands, any local or regional motions should be 
determined. Through a joint effort of the University of Hawaii's 
Institute of Geophysics and Institute for Astronomy, such 
motions are monitored through repeated geodetic laser surveys 
between the LURE Observatory and selected points on Maui and 
neighboring islands, repeated first order levelings of key 
points on Maui, and repeated gravimetric surveys on Maui and 
Oahu. Also used are ocean tide gauges, seismic monitoring of 
crustal activity, and tilt meter monitoring of changes in the 
local vertical (Carter et al,, 1977). 

This report describes the results of gravimetric measure- 
ments made in 1976-1978, compares the values with older measure- 
ments, and analyzes the reliability of the results. If the 
above procedure is done yearly, any significant changes in 
absolute gravity values could be the signal for further inves- 
tigations. Utilization of periodic first order leveling reveals 
whether changes in the absolute gravity value of a station are 
caused by a change in elevation or a change in the distribution 
of masses beneath the station. In the case of no measurable 
change in elevation it would be assumed that the distribution 
of masses had changed, and other geophysical investigations 
would be called for. Since yearly first order leveling would 
be costly, the gravity survey can also be used as a rough 
leveling survey to interpolate the periodic first order 
leveling results. 
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II. INSTRUMENTATION 

Gravity measurements were made with either of two gravity 
meters. The first was a La Coste and Romberg Model G Geodetic 
Gravity Meter No. 1; later measurements were made with a La 
Coste and Romberg Model G Geodetic Gravity Meter No. 144. The 
two gravity meters are similar in construction. Both have a 
reading precision of ± 0.01 milligal. The meters have a weight 
on the end of a horizontal beam supported by a zero-length 
spring. Their measuring system uses a measuring screw and 
lever system. Both meters have an internal heating unit to 
keep the temperature of the measuring system constant at 52.2°C 
for the G-l and 5l.85*C for the G-144. 

Operating procedure is the same for both meters. After 
the instrument is leveled, the measuring screw ic turned to 
bring the crosshair viewed in the telescope eyepiece to the 
same spot on the scale for each reading. This means that for 
each measurement the weight is nulled at the same point. The 
reading taken is then a measure of the change of length of 
the spring between gravity points, and, therefore, of the 
change in gravity acting upon the weight. The correct calibra- 
tion factors must be applied to the dial reading obtained at 
each point, before earth tide corrections and instrument drift 
correction are applied. After these corrections have been 
made differential gravity values in milligals can be found 
between each measured point. 

The theoretical earth tide corrections were computed using 
a computer program assembled by H. C. Marsh of the Hawaii 
Institute of Geophysics in 1973. The program utilized the 
Longman equations written for computers. The vertical component 
of the lunar tidal force is given as (Woollard et al., 1973); 


GM » (3 cos 2 _ - 1) + (5 cos 3 9 • 3 cos g) (1) 

d d 

where p, is Newton's gravitational constant, M is the mass of 
the moon, r is the distance from the observation site to the 
center of the earth, d is the distance from the center of the 
earth to the center of the moon, and 9 is the zenith angle of 
the moon. The vertical component of the solar tidal force is 
given as ; 


GS * (cos 2 0-1) 

D 


( 2 ) 


where S is the mass of the sun, D is 
center of the earth to the center of 
angle of the sun, and ^ and r are as 
component is, therefore, given as: 


the distance from the 
the sun, 0 is the zenith 
above. The total vertical 
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GO • GM + GS (3) 

To investigate the drift pattern of the La Coste and 
Romberg Model G Geodetic Gravity Meter No. 1 two sets of tests 
were made. During the first test period the instrument was 
placed on the floor of Room 109A In the Hawaii Institute of 
Geophysics building and was carefully leveled. Before the 
start of the test the instrument had been kept continuously 
on heat for about one week. The test period lasted 174 hours 
and five to seven readings were made daily after careful 
leveling. The clamp of the horizontal beam was kept open 
continuously during the entire test period, the room tempera- 
ture was kept constant at 23°C during the day and at about 
25°C at night, and the temperature of the measuring system 
was a constant 52.2 a C. In Table la, dates, time and number 
of observations are given, together with the daily mean and 
the maximum daily deviation from the mean. 

Twenty-one days after the end of the first set of tests 
the instrument was placed on the same location as before and 
identical testing procedure was repeated lasting this time 
295 hours. The results are given in Table lb. Between the 
two sets of tests the gravity meter was transported to Maui 
via air, used for seven days and transported around the island 
by car on very rough roads, and finally returned to Oahu via 
air. In spite of the rather small combined travel and non 
travel drifts shown in Tables la and lb we decided to close 
every survey loop to its initial point within one day. As an 
additional precaution, before opening the clamp of the horizon- 
tal beam after the gravity meter was properly leveled at gravity 
points, we turned the measuring screw to a precomputed or 
estimated correct position to prevent the crosshair bouncing 
within its mechanically limited tolerance. 
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TABLE la. Drift Pattarn of G-l Gravity Metar 


Data of 
Obaarvatlona 

Time 

No. of 

Observations 

Dally Mean, 
mgal 

Maximum 

Dally 

Difference , 
meal 

7.22-76 

0900-1600 

7 

2132.016 

+ 0.014 

7.23-76 

0900-1300 

5 

2.032 

- 0.012 

7.26-76 

1000-1630 

5 

2.094 

- 0.024 

7.27-76 

1200-1600 

5 

2.056 

- 0.005 

7.28-76 

1100-1545 

5 

2.064 

+ 0.016 

7.29-76 

1000-1500 

5 

2.068 

+ 0.012 


TABLE 

lb . Drift 

Pattern of G-l 

Gravity Meter 


Date of 
Observations 

Time 

No . of 

Observations 

Daily Mean, 
mgal 

Maximum 

Daily 

Difference , 
mgal 

8.19-76 

0900-1700 

9 

2131.890 

+ 0.010 

8.23-76 

1000- 1500 

5 

1,992 

+ 0.008 

8.24-76 

1000- 1600 

6 

1.953 

- 0.013 

8.25-76 

0800-1600 

7 

1.944 

+ 0.016 

8.26-76 

1100- 1500 

4 

1.940 

- 0.010 

8.30-76 

0900-1600 

7 

1.973 

+ 0.017 

8.31-76 

0900- 1600 

5 

1.980 

- 0.010 
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III. GRAVITY NETWORKS 

The network on Maui conalats of nine gravity points aelectad 
so that it includea five stations along the gravity calibration 
line from Kahului Airport to the summit of Mt. Haleakala. These 
are: Station 1 (Kahului Airport secondary base station), station 

3, station 5 t station 15 (Haleakala National Park Headquarters), 
station 21 (Meea Observatory). The calibration line was estab- 
lished on Maui in 1964-1965 by the personnel of the Hawaii 
Institute of Geophysics to be used in testing the calibration 
constant of gravity meters. In addition, observations were 
made at new pointB in Hana Bay, La Perouse Bay, and Pier 2, all 
of them located next to tidal gauges. An extra gravity point 
was added to Hana Airport to be compared with surveys made at the 
same site in 1961 by the 1381st Geodetic Survey Squadron, Air 
Photographic and Charting Service, United States Air Force. 
Location of these gravity points can be st»en in Fig. 1, and 
further descriptions follow. 

1. Gravity station 1 at Kahului Airport is in the terminal 
building along the wall near the airport manager's office. The 
station is marked by a small brass disk on the floor. 

2. Pier 2 station is at the western end of Pier 2 next to 
the USCGS bench mark #8 near the Kahului tidal gauge. It is 
marked by a Phillip's head screw In a lead insert drilled into 
the concrete. 

3. Gravity station 3, located on the north side of Highway 
37 approximately 0.3 miles southeast of milepost 3, is marked 

by a 12-inch-square stepping-stone set in concrete at ground 
level . 


4. Gravity station 5 at the junction of Highways 37 and 
377 is in the middle of a triangular section of ground outlined 
by the junction of the two highways, and is marked by a 12-inch- 
square concrete block set at ground level. 

5. Station 15 at the Haleakala National Park Headquarters 
is on a stone curb approximately 20 m southwest of the entrance 
to the Park Headquarters building. Station is marked by an X 
chiseled into the stone curb. 

6. Station 21 in Mees Solar Observatory is located in the 
generator room of the solar observatory. Station is marked by 
a brass disk on the floor. 

7. Hana Bay station is on a concrete wall that leads toward 
the pier at Hana Bay. It is marked by a Phillip's head screw 

in a lead insert drilled into the concrete. The station coincides 
with the Hana Bay laser terminal, and is near the tidal gauge. 



7 


8. Ban* Airport station is located at the southaast cornar 
of an opan lanal on tha vast aida of tha taralnal building. 

9. La Parouaa Bay station, at tha and of Kihel Road by 
La Parouaa Bay, is tha southernmost of tha thraa USCGS bench- 
aarks locatad naar tha shoralina, and is naar tha tidal gauga. 

On Oahu, a nine-station network for gravity aaasuraaants 
was sat up. All of tha points vara in tha Honolulu araa, froa 
Hlckaa AFB to St. Louis High School. 

1. Tha Hlckaa AFB station is inside tha Military Air 
Coaaand tarainal building on tha floor next to tha doors leading 
froa tha U. S. Custom* inspection area. 

2. Tha intar-island terminal station is locatad at the 
Honolulu International Airport beside a stone pillar at tha 
inter-island taralnal building. The point Is aarked by a brass 
disk on tha sidewalk. 

3. Station UH 325 is on tha southarnaost part of Sand 
Island on a concrete pier at the shoreline. 

4. Station UH 324 is south of the road leading to the 
State park on Sand Island. The point is on a concrete slab 
near a wire fence approximately 20 a froa the park entrance. 

5. The Bishop Museum station, in Rooa 2 of the Bishop 
Museua annex, is aarked by a brass disk on the floor. 

6. The Office of Weights and Measures station is in the 
Hawaii State Department of Agriculture Weights and Measures 
building at 1428 S. King Street. The point is on the floor 
beside a large generator in the east rooa of the building. 

7. The Hawaii Institute of Geophysics station, outside 
H1G Rooa 108 near the wall, is aarked by a brass disk on the 
sidewalk . 


8. Station UH 47 is located on the north side of Lowrey 
Street directly across froa a church. The point is on a 12-inch- 
square concrete block next to the curb at ground level. 

9. Station UH 171, on the west side of the road leading 
into St. Louis High School, is on a concrete block next to the 
curb at ground level approximately 50 a froa Waialae Avenue. 
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I?. GRAVITY NETWORK ADJUSTMENTS 

The procedure utilised le the aaaa for each gravity net- 
vork. Aftar measurements batwaan tha gravity stations hava 
been aada and corractad for aarth tides and instrument drift, 
a valua for tha gravity diffarance batwaan atationa can ba 
calculated. Thaaa differencea can ba uaed to fora loop condl- 
tiona for tha network, which can than be put into a least 
squares conditional adjuataent process. Aftar tha adjustment 
process, a discrepancy-free network is obtained along with tha 
aost probable values of tha gravity differences between stations. 
If one of tha station is assigned as a base along with an 
associated observed gravity value, the observed gravity value 
of each point can than be calculated. 

Aftar all measured lines have been tabulated, the condition 
equations for the adjustment can be found. Table 2 gives each 
measured line, the date the measurement was made, the gravity 
difference between the terminal points i> milllgals, and the 
instrument used for the Maui gravity n#'*>tk (Schenck, 1978). 
Table 3 gives the same data for the Oal u »r.avity network 
(Schenck, 1978). Figs. 2 and 3 show the Haul and Oahu gravity 
networks, respectively, along with the number of measurements 
of each line. 

When developing the condition equations for the adjustment, 
we found it easiest to vork with only triangular loops. After 
all measured lines have been used as one side of a loop, the 
number of conditions can be found easily by looking at t^e 
number of extra times each of the lines of the loop has been 
measured. It is aost important to remember that only lines 
that have actually been measured can be used when setting up 
the condition equations. For example, the loopistation 3 - 
station 21 - Hana Bay - station 3 was observed in the following 
way (Fig. 2). Station 3 was used as the base for the loop; the 
gravity meter was set up and readings were made at that station. 
Then the gravity meter was transported by car all the way to 
the summit of Mt . Haleakala to station 21 where it was set up 
and readings made. Since there are no direct roads between the 
top of Haleakela and Hana Bay, the gravity meter was transported 
back to station 3, bypassing it and proceeding along the Hana 
Highway to Hana Bay where readings were made. The gravity meter 
was then transported back to station 3 where readings were m-de 
and the loop was closed. In this manner the gravity difference 
between station 21 and Hana Bay became independently observed. 

Another important fact to remember is that each of the 
conditions must contain at least one measurement that does not 
occur in any of the other conditions. From the 43 measured 
lines of the Maul network, there were 33 condition equations that 
eould be formed. The Oahu network contained 40 measured lines 
and 32 conditions. 



1-3 

+ 

27.45 

G-l 

1-3 

+ 

27.44 

G-l 

1-3 

+ 

27.45 

G-l 

1-3 

+ 

27.46 

G-l 

1-3 

+ 

27.45 

G-144 

1-3 

+ 

27.44 

G-l 

1-3 

+ 

27.44 

G-144 

3-5 

+ 

68.52 

G-l 

3-5 

+ 

68.53 

G-l 

3-5 

+ 

68.53 

G-l 

3-5 

+ 

68.54 

G-l 

3-5 

+ 

68.52 

G-l 

3-5 

+ 

68.55 

G-144 

5-15 

+ 

321.92 

G-144 

5-15 

+ 

321.89 

G-l 

5-15 

+ 

321.89 

G-144 

15-21 

+ 

240.66 

G-144 

15-21 

+ 

240.65 

G-l 

15-21 

+ 

240.68 

G-144 

3-H.B. 
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G-144 

i.-H.A.P. 
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3-L.P. 

- 
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105.98 
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G-144 
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+ 

663.70 
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3-15 

+ 

390.41 
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TABLE 2. Maul Gravity Lines (Cent.) 



5-78 2-L.P. - 4.84 G-144 35 

5-78 2-L.P. - 4.83 G-144 36 

5-78 21-L.P. - 668.58 G-144 37 

5-78 2-5 + 101.15 G-l 38 

5-78 1-2 - 5.16 G-144 39 

5-78 3-21 4 631.08 G-144 40 

5-78 3-21 4 631.10 G-144 41 

5-78 21-H.B. - 710.00 G-144 42 

5-78 21-H.A.P. - 700.09 G-144 43 
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TABLE 3. Oahu Gravity Lines 


Date 

Line 

Gravity Dlff. 

Instrument 

Line 1 

11-77 

HIG-II 

+ 

25.80 

G-l 

1 

11-77 

HIG-II 

+ 

25.79 

G-l 

2 

12-77 

HIG-II 

+ 

25.83 

G-l 

3 

12-77 

HIG-BM 

+ 

5.87 

G-l 

4 

11-77 

HIG-W/M 

+ 

5.22 

G-l 

5 

12-77 

HIG-W/M 

+ 

5.25 

G-l 

6 

12-77 

HIG-#47 

- 

7.89 

G-l 

7 

12-77 

HIG-#171 

+ 

3.48 

G-l 

8 

12-77 

HIG-#324 

+ 

20.04 

G-l 

9 

12-77 

HIG-I325 

+ 

27.66 

G-l 

10 

12-77 

#47-#171 

+ 

11.36 

G-l 

11 

5-78 

#47— #171 

+ 

11.40 

G-144 

12 

12-29 

147-W/M 

+ 

13.12 

G-l 

13 

5-78 

#47-W/M 

+ 

13.12 

G-144 

14 

2-78 

#47-BM 

+ 

13.74 

G-l 

15 

12-77 

#171-W/M 

+ 

1.73 

G-l 

16 

12-77 

#1 71— # 325 

+ 

24.14 

G-l 

17 

12-77 

#171-#324 

+ 

16.54 

G-l 

18 

12-77 

#171-0324 

+ 

16.51 

G-l 

19 

12-77 

#325-BM 

_ 

21.74 

G-l 

20 

12-77 

#325-11 

- 

1.85 

G-l 

21 

12-77 

#325-Hick. 

- 

0.43 

G-l 

22 

12-77 

#325-0324 

- 

7.60 

G-l 

23 

12-77 

BM-W/M 

_ 

0.63 

G-l 

24 

12-77 

BM-II 

+ 

19.92 

G-l 

25 

11-77 

II-W/M 

_ 

20.57 

G-l 

26 

12-77 

II-Hick. 

+ 

1.40 

G-l 

27 

5-78 

II-Hick. 

+ 

1.37 

G-144 

28 

5-78 

II-Hick. 

+ 

1.37 

G-144 

29 

2-78 

BM-Hick. 

+ 

21.32 

G-l 

30 

2-78 

0171-Hick. 

+ 

23.75 

G-l 

31 

GO 

l 

0171-bM 

+ 

2.38 

G-144 

32 

5-78 

BM-0324 

+ 

14.14 

G-144 

33 

5-78 

#324-1 I 

+ 

5.75 

G-144 

34 

5-78 

Hick.-W/M 

_ 

21.97 

G-144 

35 
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TABLE 3. 

Oahu Gravity Lines 

(Cont . ) 


Date 

Line 

Gravity Dlff. 

Instrument 

Line No. 

5-78 

#171-11 

+ 22.27 

G-144 

36 

5-78 

#324-Hick. 

+ 7.12 

G-144 

37 

5-78 

Hick.-#47 

- 35.06 

G-144 

38 

5-78 

#47-11 

+ 33.69 

G-144 

39 

5-78 

Hick.-HIG 

- 27.21 

G-144 

40 



Maui gravity network (from Schenck, 1978). 


Mickom AFB 



Figure 3. Oahu gravity network (from Schenck, 1978) 
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The adjustment procedure 
tial leveling network, except 
differences are equal to one. 
obtained as 


is similar to that of a differen 
that the weights of all gravity 
The number of conditions, r, is 


r * n - u ( 4 ) 

where n is the total number of observed gravity differences 
and u is the number of gravity differences needed to determine 
observed gravity values to all stations with reference to the 
selected base station. From Figs. 2 and 3 it is easy to see 
that in both gravity networks, eight gravity differences were 
needed for that determination. This is a good check to the 
method described above in determining the number of conditions 
in complex networks: On Maui, r ■ 43 - 8 * 35 and on Oahu, 

r - 40 - 8 - 32. 


From the loop misclosures the condition equations are 
f ormed : 


/ 

La 


a v 

j j 


+ w. 




0 

0 


(5) 


u 


r j v j 


+ w 


r 


0, j = 1,2,3, 


n 


where a^,bj ... are known constants (in this case +1, - 1, or 
0), v^,W 2 ... are misclosures, and v 1# V 2 ... are corrections 
that when applied to the observed gravity differences will make 
the misclosures zero. 


In Eq. 5, n is larger than r and the so-called conditional 


minimum 

must be 
in Eq . 

1 » 2 , 3 , . 


adjustment is applied. In this case £ v. = minimum 

J-l J 

satisfied simultaneously with the condition equations 
5. If these equations are multiplied by -2K^, i ■ 

..,r it is obvious that the following expression is valid: 


- 2K 


rT 

. L 


Vi 


- Wj] - 2K 2 [T 


Vj 


- W, 


• • • 


- 2K 




T 1 

) 

u 


r j v j 



minimum 


( 6 ) 
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Eq, 6 is minimized by writing the partial derivatives with 
respect to each of the residuals v equal to zero. These resid- 
uals are then solved as functions of the correlates as 
follows: 


V 1 " K l a l + K 2 b l + ••• + K r r i 


V 2 " K l a 2 + K 2 b 2 + *** + K r r 2 


+ K.a + K 0 b n + ... + K r 
n l Ti n r n 


(7) 


Substituting Eqs. 7 into Eqs. 5 and combining coefficients 
the following norma 1 equations are obtained: 


K 


V 


1 1 S l a J + K 2 Z a j b J + "• + K r i a J r J + W 1 


K. Y a b + K, Y b b + ... + K ; b r + w„ 

r L-> j j 2 


( 8 ) 


• ee e •• eee 

K i Y a r . + K. Y b r + ... + K Y r.r + w *0 
1 i- j j 2 Lj j j t L, j J r 

If the misclosures in the normal equations are taken to the 
right side of the equal sign, the normal equations then form 
an r x r square matrix in which the K's are the unknowns. The 
correlates K are found by solving the normal equations by a 
suitable method and the corrections v are then obtained from 
Eq. 7. These corrections give the observed gravity differences 
the most probable values, and make the networks discrepancy 
f ree . 


The quality of the gravity survey can be estimated as a 
byproduct of the adjustment process. The standard error of 
unit weight (standard error of one observed gravity difference) 
is obtained from either of the two following formulas: 
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[Jv^/r] 172 ; j - 1,2,3, ...n 

(9) 

[I Vi'-]"*' 1 * l < 2 > 3 >--- r 

(10) 


Since each of the networks contained a large number of 
conditions a should give a good indication of the accuracy of 
the gravity meters, of the repeatability of the observations 
of the gravity differences, and of the strength of the entire 
networks. On Maui the adjustment yielded a “ ± 0.017 mgal and 
on Oahu q m ± 0*021 mgal. 

The Fortran program GRAVNET in Appendix A was used to 
carry out the conditional adjustment described above. Instruc- 
tions for the input to the program are given in the program 
itself. The program solves the normal equations by the use 
of a Gaussian elimination process. After the adjusted gravity 
differences are computed, the program goes back and checks to 
see that the original condition equation residuals have gone 
to zero. Sample outputs of the program for the adjustment of 
the Maul and Oahu networks are in Appendix B. 

All observed gravity values on Maul and on Oahu prior to 
1964 were referenced to the International Gravity Base Network 
(tfoollard and Rose, 1963). The base station of the Hawaii 
gravity network is located at the Military Air Transport System 
(MATS) terminal at Hickam Air Force Base, Honolulu. The Potsdam 
value of this gravity base station is 978,933.70 mgal. All 
observed gravity values given in this report are referenced to 
the International Gravity Standardization Net 1971 (IGSN 71) 
(Morelli e_t al., 1974). The IGSN 71 values are obtained by 
subtracting T7.S7 mgal from the Potsdam values at the Hickam 
gravity base station. Because the terminal building had been 
renovated, the original base station could not be occupied 
during the 1976-1978 gravity survey. Therefore, another 
station in the vicinity was occupied to be included in future 
networks, and for the annual comparison of inter-island gravity 
surveys the secondary gravity base station at the Kahului 
Airport, Maui was selected as the reference station with the 
IGSN 71 value of 978,874.90 mgal. University of Hawaii gravity 
stations 171, 324, and 325 also were offset from the original 
locations and were surveyed only for the future comparisons. 

Wooliard (Woollard, 1978) has presented the change in 
theoretical gravity values brought about by the recent adoption 
of the new reference ellipsoid (Geodetic Reference System, 

1967) called GRS 67. In addition, he has evaluated the relia- 
bility of the gravity standard and Potsdam datum correction 
incorporated in the values of IGSN 71. 
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Table 4 gives all gravity stations of both networks, the 
adjusted gravity differences from the base station at Kahului 
Airport, the observed gravity values of each station, and the 
station elevation. 
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TABLE 4. Gravity Differences from Base Station 
and Observed Gravity Values 


Gravity Station 

Difference from 
Base Station, 
atal 

Observed 

Gravity, 

aa_sl 

Station 

Elevation, 

m 

1 (Kahulul Airport Base) 

0 

978,874.90 

12 

2 (Pier 2) 

5.17 

880.07 

2 

3 

- 27.43 

847.47 

146 

5 

- 95.98 

778.92 

529 

15 (Park Headquarters) 

- 417.88 

457.02 

2144 

21 (Mees Observatory) 

- 658.54 

216.36 

3042 

La Perouse Bay 

10.00 

884.90 

2 

Hana Airport 

41.54 

916.44 

18 

Hans Bay 

51.48 

926.38 

3 

Inter-Island Terminal, Oahu 

43.52 

918.42 

3 

Hlckam Air Force Base 

42.13 

917.03 

6 

UH 324 

49.29 

924.19 

4 

UH 325 

41.68 

916.58 

1 

Bishop Museum 

63.44 

938.34 

24 

Office of Weights and Measures 

64.09 

938.34 

9 

HIG 

69.32 

944.22 

22 

UH 47 

77.20 

952.10 

59 

UH 171 

65.82 

940.72 

27 
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V. COMPARISON BETWEEN THE PRESENT 
AND EARLIER GRAVITY SURVEYS 

The eastern Pacific gravity survey was accomplished in 
1961 by the 1381st Geodetic Survey Squadron (Missile) , Air 
Photographic and Charting Service, United States Air Force, 
Orlando Air Force Base, Florida, The survey was made in support 
of the Hawaiian HIRAN Project, AF 61-2, and Included base station 
surveys and observations along the level lines. All final data 
were referred to the MATS terminal base station located at 
Hickam AFB , Honolulu (USAF Report, 1963). 

Two Worden Master gravity meters were used. These instru- 
ments, number 615 and 617, are of the unstable type and have a 
three-piece quartz spring element. The effects of temperature 
on instrument drift variability and rate were minimized by means 
of the insulating vacuum flask protecting the instrument, the 
temperature variation compensation system, and the low-powered 
temperature stabilization system that provides actual temperature 
control. The quartz spring element is sealed in partial vacuum 
for protection against contamination and to minimize the effects 
of changing atmosphere. 

The “Tidal Gravity Effect Tables'* published by the Houston 
Technical Laboratories were used to determine and remove the 
effect of the attraction of the sun and moon on the gravity 
measurements. As additional data from which the quality of the 
survey could be estimated linear regressions of 12 observations 
by instrument 615 and of 13 observations by instrument 617 were 
performed at Hickam AFB during the period 6 July to 7 August, 

1961. For instrument 615 a drift rate of 0.5307 mgal per day 
was determined, together with the standard error jf unit weight 
- ± 0.11 mgal. Corresponding figures for instrument 617 were 
0.3360 mgal per day and 02 “ * 0.09 mgal. In the final computa- 
tions the mean of the observations made by the two instruments 
was used; estimated standard deviation per gravity difference 
was Ofl • ± 0.07 mgal. Three stations surveyed in 1961 were 
resurveyed in 1976-1978. Their locations and the 1961 observed 
gravity values in IGSN 71 system are: 

Kahului Airport, Maui 978,874.90 mgal 

Hana Airport, Maui 978,916.52 mgal 

Bishop Museum, Oahu 978,938.38 mgal 
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The gravity survey on Oahu resulting in the computations 
of Bouguer anomaly maps mas completed by the personnel of the 
Hawaii Institute of Geophysics in 1964* The work was done by 
using La Coste and Romberg Geodetic Gravity Meters with an 
estimated standard deviation of » * 0.02 mgal per gravity 
difference. Five stations surveyed in 1964 were resurveyed 
in 1976-1978. Their location and the 1964 observed gravity 
values in IGSH 71 system are: 


Location 

Mgal 

Inter-island terminal 

978,918.42 

Bishop Museum 

978,938.35 

Office of Weights and 
Measures 

978,939.01 

Hawaii Institute of 
Geophysics 

978,944.30 

University of Hawaii 
gravity station 47 

978,952.12 


In 196$ a calibration line was established on Maul to test 
the calibration constants of gravity meters. Originally this 
line consisted of 21 stations starting at the Kahulul Airport 
and ending at the generator room in the Mees Solar Observatory. 
From 1965 to 1978 most of the stations were destroyed, either 
by vandalism or by new road construction, and could not be 
relocated. Of the stations available, five were selected to be 
included in the Maui gravity adjustment network. The 1965 values 
were obtained as mean values of six observations made by using 
La Coste and Romberg Geodetic Gravity Meters during 1965. The 
gravity meters were: G-12, G-19, G-64, and G-65. In Table 5 

the mean gravity values and their standard deviations for the 
five stations are gi’»en along with the specific numbers of the 
gravity meters and dates of observations. 

In the 1961 and 1964-1965 gravity surveys, total instrument 
travel drift was taken to be the loop closure obtained after 
correction of observations for non-travel drift and tidal 
effects. Travel drift corrections for each observation were 
then obtained by prorating the above loop closures on the basis 
of elapsed travel time of observation from the loop base station. 


TABLE 5. Gravity Survey on Maui In 1965. Values are in *gal 
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Table 6 summarizes the observed gravity values in 1961, 
1964-1965, and 1976-1978 surveys, together with their standard 
deviations, differences between the old and present gravity 
values, and standard deviations of the differences obtained 
from 00 • ± ( 0 q 2 + o^ 2 )l/ 2 > it is obvious that only a difference 
larger than its standard deviation could warrant further inves- 
tigation as to its cause. As can be seen, only two of the 
gravity stations can be considered to have shown significant 
changes. Gravity station 15 at the Park Headquarters on Maui 
and the gravity station at the Hawaii Institute of Geophysics on 
Oahu need close surveillance. Gravity differences at the gravity 
station 3 and at the Hana Airport on Maui are too close to their 
standard deviations to be considered significant. 


TABLE 6 . Comparison of Gravity Surveys Made on the Islands of Maui 
and Oahu in 1961 , 1964 - 1965 , and 1976 - 1978 . Values are in mgal. 


24 


co 

cm cn •j in 
o o o o 
♦ • • * 

cn 

o 

• 

m 

o 

a 

© © © 
• « • 

CM 

O' 

o o o o 

© 

© 

o © o 

1— 1 
1 

U*V 

M3 

O' 

•H +1 -M +1 

+1 

+t 

+i +i +i 

niicnn 
O OHO 
• • • • 

g 

• 

H 

o 

• 

CM 00 CM 

O O O 
• • * 

a h 

HI 

o o o o 

o 

© 

o o © 

1 

+ 1 1 1 


+ 

+ + + 


u 


CO 
•H f* 

a O' 


O ' 1 

rH 


V 


I 

VO 

rv 

O' 


a> 


vO 

O' 


u 

•H 

> 

o 


o 


o 


o 

+1 


o 

•H 


ao 

o 


o 

+ 


o 

+ 


CMCMCMCMCMCMCMCM 

OOOOOOOO 


N N N N N 

o o o o o 


OOOOOOOO 

+1 +1 +1 +1 +1 +1 +1 +1 


o o o o o 

+1 +1 -N +1 +1 


O I s * 
O' o 


S CM CM VO O -O OO 

»ono"» n 


N<*)OlCO«J 
■o- O rH in cn 


-<r o 

f" CO 

co co 


N 00 N « < VO « 
co r- -a- cm co O' O' 


00 C*" «o vO 00 
H H N H fO 
O' O' O' O' O' 


00 

ts 

O' 


r*. 

o 


o 


o 

+1 


o 

O' 


CM 

m 


00 

cn 


r-* 

00 


VO 

rH 

O' 


00 

CO 

O' 


CM CM 

o o 


CM CM 

o o 


o o 

+1 +1 


o o 

+1 +1 


O' CM 
O' CM 


O CM 


00 

c n ■O’ 

O' O' 


CM O 

•n >J 

O' O' 




h cm n in 


CM 

CM 

CM CM CM 

► 


o O o o 


O 

O 

O © O 

0) 


» * . • 

3 

• 

• 

• if 

> 


o o o o 

43 

o 

o 

o o © 

u 

•Hi 






3 

CO 

2 

+1 +1 +1 +1 

ol 

+1 

+1 

+1 +1 +1 


SO 

OHOvM 


CM 

m 

rH O CM 

n 

O' 

m O' co cn 



cn 

O cn rH 

vo 

• 

.... 


• 

• 

• * • 

O' 

<r 

N 00 vo vO 


oo 

00 

O' -o CM 

iM 

r- 

•I r. in h 


rH 

cn 

cn m 

1 

co 

co r-. sr cm 


O' 

O' 

O' O' O' 

•O 

* 






vO 

00 






O' 

C"- 






rH 

O' 







00 

r>. 

O' 


o 

& 


< 

*W ' 

3 CM 


3 u 

JZ <u 

c2 2 


CO r-> 


rH OJ 



U fn 


CO CO 



OJ u 


c 3 



U 0 
Ih 4H 


iH CO 

a 

CO 


CO CO 


M V 

4J 


3 > 


OJ u 

43 


O’ >H ►, 


JJ J-l 

- oo 

CO 

*3 01 10 


0 

S ‘ H 

0) 

CO « PQ 

U 

TJ tfc 

3 OJ 

b 

OJ HP 

M 

C 

<U 3 

3 

PS o 0> 

O 

CQ U 

CO 

CO 

CO 

O. 

rH *H 

3 IJH 

<0 

M CO 3 

Cl 5n 

CO < 

s o 

OJ 


CM 


- - 0 ) _ 

co oj J* < 
fc S3 « 
'-'Cu co 
c 
« 
PC 


3 


m h o 

co in iH CM i-l 


co 

3 

PS 


C *H 


B nt m 

CO CM CM 

j* m cn 

u 

PS 33 
PC 33 S 3 


a a) 
o o 
jc h 

(0 t*H 
•H ’4H 

« o 


33 

P3 



25 


REFERENCES 

Carter, W. E., E. Berg, and S. Laurila. The University of 
Hawaii Lunar Ranging Experiment Geode t ic -Ge ophys ica 1 
Support Programme. Phil. Trans. R. Soc. Lond. A., 
v. 284, p. 451-456, 1977. 

Morelli, C., C. Gantar, T. Honkasalo, R. K. MConnel, J. G. 
Tanner, B. Szabo, U. Uotila, and C. T. Whalen. The 
International Gravity Standardization Net 1971. Spec, 
Pub. 4, Inti. Assoc. Geodesy, 194 pp., Paris, 1974. 


Schenck, B. 

S tudy . 
Hawaii , 


E. High Precision Geodesy Applied to Crustal 
MS Thesis, Dept. Geology and Geophysics, Univ. 
Honolulu, 111 pp., 1978. 


United States Air Force. Report of Results, Project AF 61-2 
Eastern Pacific Survey. Vol. IV, June 15, 1963, 111 pp. 
1370th Photo-Mapping Wing, Air Photographic and Charting 
Service, USAF, Orlando Air Force Base, Florida. 

Woollard, G. P. The new gravity system: Changes in the 

international gravity base values and anomaly values. 
Geophysics, 1978 (in press). 


Woollard, G. P., H. C. Marsh, and R. L. Longfield. A Comparison 
of Observed and Predicted Earth Tide Observations between 
Alaska and Mexico and across the Zone of Crustal Transi- 
tion between Nebraska and Utah. Hawaii Instit. Geophys. 
Rept. HIG-73-20, Univ. Hawaii, Honolulu, 1973. 

Woollard, G. P., and J. C. Rose. International Gravity 

Measurements. Spec. Pub., Soc. Explor. Geophys., Tulsa 
Oklahoma, 518 pp. 1963. 


PRECEDING PAGE BLANK NOT FNJVKD 


27 


APPENDIX A 


COMPUTER PROGRAM FOR 


GRAVITY NETWORK ADJUSTMENT 
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cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c 

c **** CRAVNET ***» 

C PROGRAM TO ADJUST CR.1V I TY NETWORK BY A LEAST SQUARES METHOD » 

C WRITTEN BY BRUCE SCHENCK. JAR. 1973. 

C 

C THIS PROCRAM TAKES AS INPUT MEASURED GRAVITY DIFFERENCES IR 

C MILL I CALS BETWEEN POINTS OF A NETWORK. FROM THE NETWORK LOOP 

C CONDITIONS ARE FORMED. AND THE MEASURED DIFFERENCES ARE THEN 
C ADJUSTED BY LEAST SQUARES SO THAT EACH LOOP OF THE NETWORK HAS 

C NO CLOSURE ERROR. AS THE PROGRAM IS WRITTEN NOW 89 MEASURED 

C DIFFERENCES CAN BE INPUT, AND 50 CONDITIONS CAN BE SOLVED. HOV- 

C EVER WITH A FEW CHANCES THESE VARIABLES CAN BE INCREASED. FOR 

C EACH NETWORK TO BE ADJUSTED THE CONDITIONS IN THE PROCRAM MUST 
C BE CHANGED ACCORDINGLY, ALONG WITH THE INITIALIZATION OF THE 
C CORRELATE MATRIX. 

C INPUT TO THE PROGRAM IS AS FOLLOWS : 

C CARD *1: HEADER CARD. ANYTHING ON THIS CARD WILL BE PRINTED AT 
C THE TOP OF THE OUTPUT. 

C CARD *2: *NO“. THE NUMBER OF MEASURED DIFFERENCE TO BE READ IN 
C BY THE PROGRAM, IN AN 12 FORMAT STARTING IN COLUMN *4. 

C CARD *3: "N*, THE NUMBER OF CONDITIONS PRESENT IN THE ADJUSTMENT 
C IN AN 12 FORMAT STARTINC IN COLUMN' *3 . 

C CARD *4-0N: COLS. 1-12: TERMINAL POINTS OF MEASURED DIFFERENCES 

C IN A 2A6 FORMAT s 

C COLS. 15-21: MEASURED GRAVITY DIFFERENCE IN MILLI- 

C CALS, IN AN F7.0 FORMAT; 

C COLS. 25-28: WEIGHT OF MEASUREMENT IN AN F4.0 FORMAT 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
c*»« THIS SECTION INITIALIZES VARIABLES AND READS INPUT. 

DIMENSION ALIKE! OA) ,EX( 12) .XLINEC 160) , W(50) .VT(O0) ,R( 80.30) . 

ICC 30, 30) ,S<30> ,QC SO) .WWC30) , V(80> 

NN«0 

READ CO, 100) EX 

100 FORMAT ( 12A6) 

READ (8, 101) NO 

101 FORMAT <3X. 12) 

READ C 8, 126) N 

120 FORMAT ( 2X, 12 > 

READ (8,102) (XLINEC 2*1-1) ,XLINEC2*I) , ALINE( I) ,VT( I) . 1*1, NO) 

102 FORMAT ( 2A6 , Tt5 , F7. 0 , T25 , F4. 0) 

WRITE C9. 100) EX 

WRITE (9, 121) 

121 FORMAT </, ’ OBSERVED GRAVITY DIFFERENCES:’) 

119 CONTINUE 

WRITE (9 103) 

WRITE ( 9 ! 104) ( I , XLINEC 2* I- 1 > ,XLINE(2*I) .ALINEC I) ,WT( I) , 1*1. NO) 

IF CNN . EQ. 1) GO TO 122 

103 FORMAT ( ’ ,T6, 'LINE: ' .T20, ’CRAV. DIFF. T35 WEIGHT: • ) 

104 FORMAT C M2, T6 ,2A6 , T20 ,F7. 2, T35 , F4. 1) 

DO 105 1*1.30 

105 wc i:**o. 

123 CONTINUE 

C*** THIS SECTION CONTAINS THE CONDITIONS FOR THE ADJUSTMENT. 

C****> NOTE: THE CONDITIONS MUST 3E CHANGED FOR EACH DIFFERENT NETWORK. 
C*** OAHU CONDITIONS PRESENT* 

WC I > —ALINEC 25) + ALINEC 33) ♦ALINEC 34) 

WC2) *-ALINE( 25) -ALINEC 32) ♦ALINEC 36) 

WC3) *- ALINEC 25) -ALINEC 15) ♦ ALINEC 39 > 

VC 4) —ALINEC 25) +ALINEC 24) -.ALINEC 26) 

WC 5) —ALINEC 25) -ALINEC 4) ♦ ALINEC 1) 

WC 6) *- ALINEC 25) -ALINEC 4) + ALINEC 2) 


^ c 

!HTt- 


oowooooyoooooooooooououuuouo 


V< 7) —ALINE! 25) - AL TNE! 4) +A L I NEC 3) 

V< 8) —ALINE! 25) +ALINE! 30) -ALINE! 87) 

W!9)—ALINE(25>+ALINE( 30) -ALINE! 23) 

W< 10) —ALINE! 25) + ALINE! 30) -ALINE! 29) 

V! 1 1) —ALINE! 25) -ALINE ( 20) + ALINE! 21) 

W< 12) —ALIKE! 26) -ALINE! 13) +ALINE! 39) 

13) —ALINE! 26)-ALINE! 14) +ALINE! 39) 

V( 14) —ALINE! 26)-ALINE! 16) +ALINE! 36) 

W< 15)— ALINE!26)-ALINE!5>+ALINE( l) 

\K 16>— ALINE! 26) -ALINE! 6) +AL IKE! 2> 

W( 17) —ALINE! 26)-ALINE(35)-ALINE<27) 

W( 18) —ALINE! 39) ♦ ALINE! U>+ALINE<36) 

W< 19) — ALINE! 39) ♦ ALINE! 12) +ALINE! 36) 

W( 20) —ALINE! 39 ) -ALINE! 38) -ALINE! 28) 

V! 21 >— ALINE! 39) -ALINE! 7) + ALINE! 3) 

W< 22) —ALINE! 1 ) + ALINE! 8)+ALINE! 36) 

W( 23) —ALINE! 2) + ALINE! 9) + ALINE! 34) 

W! 24) — ALINE! 3) +ALINE! 10) ♦ ALINE! 21 ) 

W( 25) —ALINE! D-ALIKE! 40)-ALINE!29) 

W! 26 )» -ALIKE! 36)+ ALIKE! 17) +ALINE! 2 1 ) 

W! 27) —ALINE! 36) +ALIKE! 18) + ALIKE! 34) 

W< 28) »- ALIKE! 36) + ALIKE! 19 ) +ALIKE! 34) 

V! 29) --ALIKE! 36) +ALIKE! 3 1) -ALIKE! 28) 

W! 30) « - AL I KE! 34 ) -AL I KE! 23 ) + ALIKE! 2 1 ) 

W( 3 l ) * - AL I KE( 34 ) +AL IKE! 37) -ALINE! 29 ) 

W! 32) --ALIKE! 21) + ALIKE! 22) -ALIKE! 27) 

C *** MAUI CONDITIONS PRESENT! 

IF !NN .EQ. 1) CO TO 124 

C*** THIS SECTION COMPUTES THE CONDITION RESIDUALS. 

WRITE (9 l 06) 

106 FORMAT <✓,’ CONDITION RESIDUALS BEFORE ADJUSTMENT: * ) 

WRITE ! 9 , 107) I I,V( I) , I«1,N) 

107 FORMAT I 13, * : ' .F5.2) 

DO 150 I-l.N 

150 VW! I)*W( I) 

C*** THIS SECTION INITIALIZES THE CORRELATE MATRIX. 

DO 108 1-1.50 
DO 108 J-1,80 

108 R! J. I) *0. 

C-»-«> NOTE: THE INITIALIZATION Ofc THE CORRELATE MATRIX MUST BE CHANCED 
C FOR EACH NETWORK. 

C»*» CONDITION * 1 j INITIALIZE 
R!25. 1> — 1. 

R! 33. 1 ) *♦ 1 . 

R! 34, 1 ) » + 1 . 

C*** CONDITION *2; INITIALIZE 
R( 25.2)--l. 

R! 32.2) — 1. 

R! 36 , 2) — l . 

C*** CONDITION *3i INITIALIZE 
R( 25 , 3 ) * - 1 . 

R( 18,3) — 1. 

R! 39 , 3) ■ + 1 . 

C*** CONDITION #4; INITIALIZE 
R! 25 , 4) *- l . 

R< 24, 4) « + l . 

R! 26 , 4) » - 1 . 

C*** CONDITION *5; INITIALIZE 
R! 25 , 5 ) - - 1 . 

R! 4 . 5 ) — 1 . 

R! 1 . 5) * + 1 . 

C*** CONDITION *6 ; INITIALIZE 
R( 25 , 6) s - 1 . 

R! 4 . 6 ) - - 1 . 
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R<2,6)«+1. 

C*** CONDITION *7; INITIALIZE 
R(25,7)*-l. 

RC4.7)*-!. 

R< 3.7) «+ 1 . 

C*** CONDITION *8; INITIALIZE 
R( 25 , 8) »- 1 . 

R< 30, 8) ** l . 

R<27,8>*-1. 

C*** CONDITION *9| INITIALIZE 
R(25,9)*-i. 

R(30,9)»+l. 

R(28 9)s—i, 

C*** CONDITION *10; INITIALIZE 
R(28, 10) •-!. 

R(30, 10)*+!. 

C*** CONDITION *11: INITIALIZE 
R< 25 , 11 ) « - 1 . 

R( 20, 1 1) t . 

RC21, 1 l)»+l. 

C*»* CORD IT! OR *12; INITIALIZE 
R<26, 12)=-M. 

R( 13, I2)«-l. 

R( 39 , 12) *+■ 1 . 

0*sb CONDITION *13: INITIALIZE 
R( 26 , 18) *+1, 

R( 14, JS)*-1. 

R( 39 13) * -t*l . 

C*** CONDITION *14: INITIALIZE 
R(26, 14) a + 1 . 

R( 16, 14) *- 1 . 

R( 36, 14) »+l. 

C*** CONDITION *15; INITIALIZE 
R( 26 , 15) » + l . 

R( 3. 15) *- 1 . 

R( 1. 15) =* I . 

C*** CONDITION *16: INITIALIZE 
R( 26 , 16) •* I . 

R( 6 . 16)*-1 . 

R<2. 16) =+■ l , 

C*** CONDITION *17; INITIALIZE 
A( 26 , 17) *+I. 

R( 33 , 17) *- l . 

R( 27, 17) »- 1 . 

C**» CONDITION *18; INITIALIZE 
R(39, 18) *- 1 . 

R( II, 18)« + i. 

R( 36 . 18) *+ 1 . 

C*«* CONDITION *19, INITIALIZE 
R( 39 , 19) ■- 1 . 

R( 12. 19)» + 1. 

R( 36 , 19)**1. 

C*** CONDITION *29; INITIALIZE 
R( 89,20) «-l . 

R( 33,20) *-l . 

R(28,20)»-l. 

C*** CONDITION *21; INITIALIZE 
R( 39,21) *-l . 

R(7,21)*-1. 

R( 3, 21 ) *+ l . 

C*** CONDITION *22; INITIALIZE 
R( 1.22) *-» I . 

R(8,22)=+l. 
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R< 34,22)*+!. 

C»tt» CONDITION *23 t INITIALIZE 
R<2.23>«-1. 

R<9.23)*+1. 

R< 34,23)*+!. 

Cwmi CONDITION *24 t INITIALIZE 
R(3.24>*-1. 

R( 10,24) ■+ 1 . 

R(2l ,24) *+I. 

C*** CONDITION *23; INITIALIZE 
RU.2S) — 1. 

R( 40,23)*-!. 

R(29 23)*-l. 

C*** CONDITION *26 t INITIALIZE 
R(36,26) ■-!. 

R< 17,26) *+l. 

R( 21 ,26) * + 1 . 

C**» CONDITION *27* INITIALIZE 
R(36 .27) »- 1 . 

R< 18,27) *♦ 1 . 

R(34,27)*+l. 

C*w CONDITION *28* INITIALIZE 
R( 36,28) *- 1 . 

R( 19,28) *♦ l . 

R( 34,28) *+l. 

C*** CONDITION *29* INITIALIZE 
R( 36 , 29 ) * - 1 . 

R< 31 , 29) *+l . 

R( 28,29) »-l. 

C**» CONDITION *30* INITIALIZE 
RI34.30) *-l. 

R(23,30) *-l . 

R(2l ,30) . 

C*** CONDITION *31 * INITIALIZE 
R< 34,31) «-i. 

R( 37, 31 ) *♦ 1 . 

R< 29 , 3 1 ) ■ - 1 . 

C**w CONDITION *32* INITIALIZE 
R(2l,32)*-1. 

RC22,32)*+1. 

RC 27 32) *• l • 

C*** THIS* SECTION SETS UP THE NORMAL EQUATIONS. 

DO 109 1*1,50 
DO 109 J* I , 50 

109 C( l, J>*0. 

DO 110 K«l,N 

DO 110 J* 1 ,N 

DO 110 1*1, NO 

S( I> *R( I , K)*R( I.Jl/'VTC I) 

110 C( K, J) *C( K. J)*S( I) 

C*** THIS SECTION CONSISTS OF A GAUSSIAN ELIMINATION PROCESS TO SOLVE 
C*** THE NORMAL EQUATIONS FOR THE CORRELATES. 

IF (N ,NE. 1) CO TO 4 
IF (C( 1, 1> .EQ. 0.) GO TO 3 
Q( 1)*W( D/CC 1, 1)«(-1.) 

GO TO 202 

3 GO TO 203 

4 NLESS 1 * N- 1 

DO 13 I* 1, NLESS 1 
BIG* ABS( C( I. I>) 

L* I 

IPLUSl* 1+ 1 
DO 6 J* IPLUSl, N 

IF ( ABS( CC J, I) ) . LE. BIG) GO TO 6 
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6 continue 

IF <B!C ,5E. 0.) CO TO a 
GO TO 203 

8 IF (L .EO. I) CO TO II 
DO 10 J«l,5 
TEMP-C(L.J) 

C(L, J)*C( I , J) 

10 C( I.JMTEMP 
TEWF*W(L) 

W(L) «1f< I) 

V<I)*TEMP 

11 DO 13 J»IPLUS1,5 
OOOT»C(J. I)/CC 1. 1) 

DO 12 K«lPLUSl,N 

12 C(J.©*C< J,©-OUOT*C( I,D 

13 WC J) -W( J>-O0OT*W< I) 

IF (C(N,5> .RE. 0. > CO TO 13 
CO TO 203 

15 Q<N>*V15)/C(5,5) 

I-Pf-1 

to sun-0. 

IPIUS1*I+1 
DO IT J* IPLUS1 , IT 

17 SUM*SUM+C( I,J)*a< J) 

Q< I)*(W( I) -SUM) /Cl 1,1) 

1*1-1 

IF ( I .CT. 0.) CO TO 16 
DO 18 KK»1,5 

18 Q( KJO *-a(KK) 

202 WRITE (9,201) ( I ,Q( I) , I* 1,R) 

201 FORMAT (’ CORRELATES: ’ ,/.25( 13, * t * ,F15.3,/) ./) 

CO TO 203 

203 WRITE (9,204) 

204 FORMAT (’ ***ERROR**»**> IN CORRELATE MATRIX*) 

CO TO 999 

205 CONTINUE 
SUMW«0. 

DO (11 1*1,5 

1 1 1 SOMW»SUMW-a< I)*WV( I) 

C*** THIS SECTION COMPUTES THE CORRECTIONS TO EACH MEASUREMENT. 

DO 112 1*1,80 

112 V( I)*0. 

DO 113 1*1,50 
DO 113 J* 1 , 5 
W*R( I,J)*0(J) 

113 V( I)*V( I)*W 
DO 114 1*1,50 

114 V( I>«V( I)/WT( I) 

WRITE (9,113) ( I , V{ I) , I» 1,50) 

113 FORMAT (’ CORRECTIONS TO LINES: * 50( 13, * *• * ,F7. 3,/) ) 

C*** THIS SECTION COMPUTES SUM W , A CHECK 05 SUM W, AND MU. 

SUMW2*0. 

DO 116 1*1,50 
GI«V( I>**2 

116 SUMW2*SUMW2*GI 
AMU*SCRT( SUTIW2X5) 

WRITE (9,117) SUMW.SUMW2.AMU 

117 FORMAT ( 1 SUM W* ' .F10.3./, • CHECK 05 SUM W« * ,F18.3,/’. ' HU*',F10. 
13./) 

C*** THIS SECTION COMPUTES EACH ADJUSTED MEASUREMENT. 

DO 118 1*1.50 

118 ALI5E( I) *ALINE( I)+V( I) 
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NN»NN-M 
WRITE (9, 120) 

120 FORMAT < • CORRECTED GRAVITY DIFFERENCES! ' ) 

GO TO 119 
122 CONTINUE 

C*» THIS SECTION COMPOTES THE CONDITION RESIDUALS AFTER ADJUSTMENT. 
GO TO 123 

124 CONTINUE 

UNITE <9, 125) 

125 FORMAT </,’ CONDITION RESIDUALS AFTER ADJUSTMENT:’) 

UNITE (9,107) ( I ,U( I) , I" 1 ,N) 

909 STOP 
END 
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HAUI CRAVITY NTTVORT 
OBSERVED GRAVITY DtfTEREftCESi 


# 

LIRE: 

GRAV. Dirr.i 

WEIGHT: 

1 

*1 TO *3 

27.45 

1.8 

2 

•1 TO *3 

27.44 

1.8 

3 

*1 TO *3 

27.45 

1.8 

4 

*1 TO *3 

27.44 

1.8 

5 

*1 TO #3 

27.45 

1.8 

6 

*1 TO *3 

27.44 

1.9 

7 

*1 TO *3 

27.44 

1.9 

8 

#3 TO *5 

68.52 

1.9 

9 

*3 TO *3 

68.53 

1.9 

10 

*3 TO *5 

68.53 

1.9 

11 

*3 TO *5 

68.54 

1.9 

12 

*3 TO *3 

68.52 

1.9 

13 

*3 TO *5 

68.55 

1.9 

14 

•3 TO *15 

321.92 

1.9 

15 

*5 TO *15 

321.89 

1.9 

16 

*5 TO *15 

321.89 

1.9 

17 

*15 TO *21 

240.66 

1.9 

18 

*15 TO *21 

240.65 

1.0 

19 

*15 TO *21 

240.68 

1.9 

20 

*3 TO I.B. 

-78.92 

1.0 

21 

H.8. TO SAP 

9.98 

1.9 

22 

*3 TO HAP 

-69.01 

1.0 

23 

*3 TO HAP 

-68.99 

1.9 

24 

*3 TO L.P. 

-37.45 

1.9 

25 

*3 TO L.P. 

-37.46 

t.O 

26 

*1 TO L.P. 

-9.98 

1.0 

27 

*1 TO L.P. 

-4.97 

1.0 

28 

*5 TO L.P. 

-105.98 

1.9 

29 

*5 TO L.P. 

-105.90 

1.0 

30 

*15 TO L.P. 

-427.07 

1.0 

31 

*15 TO L.P. 

-427.91 

1.9 

32 

*2 TO *21 

663.70 

1.8 

33 

*3 TO *15 

390.41 

1.0 

34 

*2 TO *3 

32.64 

1.9 

35 

*2 TO L.P. 

-4.84 

1.0 

36 

*2 TO L.P. 

-4.83 

1.0 

37 

*2! TO L.P. 

-668.58 

1.9 

38 

*2 TO *5 

101. IS 

1.0 

39 

*1 TO *2 

-5. 16 

1.0 

40 

*3 TO *21 

63 1 . 08 

1.0 

41 

*3 TO *21 

631. 10 

1.0 

42 

*21 TO H. B. 

-719.09 

1.0 

43 

*21 TO HAP 

-700.09 

1.9 



COITD IT 10W RES I DUALS BEFORE ADJUSTMENT: 
is 9.02 


2 : 

3: 

4? 

5: 

4: 

7: 


9.03 
O.Ol 
9.03 
O.Ot 
9.01 
0.00 
8 :-0.0S 
95-9.02 
10: -0. 05 
11 .- 0.02 
12: -0. 03 
13: -0. 02 
145-0.03 


*»»&«■* 


VI 

lxc< 


it p 

n 

P45E NATIONALLY BLAf!.£ 


15 i ••.63 
141 - 6.01 
17 : *.01 
1S> 0.00 
iOi #. 02 

201-0.01 
21 1 •.•• 
22 i- 0 . 0 l 
23! -0. 04 
24) ••.OS 
25 i -• . 02 
201 - 0.01 
27t 0.00 
28: -0.04 
24: -0.05 
30-0.02 
31: -0.03 
32t -0.0* 

33: -0.05 
34! 0.00 
35.-0.02 
CORRELATES: 


n 

0.077 

2i 

-0.008 

3: 

-0.047 

4! 

-0.043 

Si 

0.023 

0> 

•0.034 

7 ! 

0.035 

8. 

0.013 

Ot 

• .•OS 

10t 

0.015 

lit 

0.005 

12! 

-•.005 

13i 

• .005 

14! 

0.415 

IS; 

• .015 

10! 

0.010 

17! 

0.000 

18: 

0.000 

14: 

-0 . 0 1 0 

20> 

0.019 

21! 

-0.020 

22; 

-0.034 

23: 

0.083 

24 ; 

0.003 

25; 

0.017 

2«: 

-0.021 

27: 

0.004 

28i 

0.004 

24: 

0.004 

30: 

0.014 

31* 

-0.014 

32: 

0.033 

33: 

0.003 

34: 

-0. 010 

35: 

0.024 


CORRECT! ORS TO LIHES* 
l! 0.005 
2: 0.015 

3: 0.005 
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4« 

-4.498 

8s 

0.445 

4s 

4.413 

?! 

4.418 

8! 

4.414 

9s 

4.444 

14s 

4.444 

Ut 

-4. 414 

12s 

4.414 

13s 

-4.024 

14s 

-4.421 

18 s 

4.449 

14 s 

9.449 

17s 

9.444 

18s 

4.414 

19s 

-4.414 

24s 

-4.411 

21 s 

-4.436 

22s 

4.923 

23s 

9.943 

24s 

-4.443 

25: 

9.447 

24s 

-4.019 

27 s 

-4.429 

28 s 

-0.003 

29s 

-0.093 

SO: 

-0.013 

31 s 

9.02? 

32: 

0.013 

33: 

0.019 

34: 

-0.020 

38: 

0.006 

34 s 

-0.044 

37: 

9.034 

38: 

0.000 

39: 

-0.005 

44: 

4.013 

41: 

-0.907 

42: 

-0.024 

43: 

0.019 




sun w* 0.019 
check on sun w* 0.010 

MU« 0.017 

CORRECTED CHAV ITT DIFFERENCES: 


m 

LINE: 


chav. Dirr.s 

WEIGHT 

1 

m 1 

TO 

'3 

27.45 

1.9 

2 

m ) 

TO 

'3 

27. 4fi 

1.9 

3 

* } 

TO 

'3 

27.45 

1.9 

4 

# 1 

TO 

'3 

27.43 

1.4 

5 

m 1 

TO 

'3 

27.45 

1.0 

6 

• | 

TO 

'3 

27.45 

1.0 

7 

m J 

TO 

*9 

27.45 

1.4 

8 

'3 

TO 


68.83 

1.0 

9 

*3 

TO 


68.53 

1.0 

10 

*3 

TO 


68.53 

1.0 

11 

'3 

TO 

•9 

68.53 

1.9 

12 

*3 

TO 


68.53 

1.0 

13 

*3 

TO 

'5 

68.53 

1.0 

14 

*9 

TO 

'15 

321.90 

1.0 

15 

*9 

TO 

'15 

321.90 

1.0 

16 

*9 

TO 

'15 

32 1 . 90 

1.0 

17 

*15 TO '21 

244.66 

1.0 
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18 

*13 TO *21 

240.66 

1.0 

14 

*13 TO **21 

240.66 

1.0 

28 

*3 TO 0 E. 

-78.93 

1.0 

21 

B.B TP BAP 

9.94 

1.0 

22 

*3 TO HAP 

-68.99 

1.0 

23 

*3 TO BAP 

-68.99 

1.0 

24 

*3 TO L.P. 

-37.43 

1.0 

23 

*3 TO L.P. 

-37. 45 

1.0 

26 

*1 TO L.P. 

-10.00 

1.0 

27 

*1 TO L.P. 

-10.00 

1.0 

28 

*5 TO L.P, 

-103.98 

1.0 

29 

*3 TO L.P. 

-103.98 

1.0 

38 

*13 TO L.P. 

-427.88 

1.0 

31 

*13 TO L.P. 

-427.88 

1.0 

32 

*2 TO *21 

663.71 

1.0 

33 

*3 TO *15 

390.43 

1.0 

34 

*2 TO *8 

32.62 

1.0 

33 

*2 TO L.P. 

-4.83 

1.0 

36 

*2 TO L.P. 

-4.83 

1.0 

37 

*21 TO L.P. 

-668. 33 

1.0 

38 

*2 TO *5 

101. 15 

1.0 

39 

*1 TO *2 

-3. 17 

1.0 

48 

*3 TO *21 

631.09 

1.0 

41 

*3 TO *21 

631.09 

1.0 

42 

*21 TO H.B. 

-710.02 

1.0 

43 

*21 TO HAP 

-700.08 

1.0 

COTO IT I ON RESIDUALS AFTER 

ADJUSTMENT: 

1 

8.08 



2 

8.00 



3 

8.88 



4 

0.80 



3 

0.80 



6 

0.00 



7 

0.00 



8 

0.80 



9 

0.00 



18 

0.00 



11 

0.00 



12 

0.00 



13 

0.00 



14 

8.00 



13 

0.00 



16 

0.00 



17 

0.00 



18 

0.00 



19 

0.00 



20 

0.00 



21 

0.00 



22 

0.00 



23 

0.00 



24 

0.00 



23 

0.00 



26 

0.00 



27 

0.00 



28 

0.00 



29 

0.00 



38 

0.00 



31 

0.00 



32 

0.00 



33 

0.00 



34 

0.00 



33 

0.00 
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OAHU GRAVITY NETWORK 
OBSERVED GRAVITY DIFFERENCES: 


• 

LINE: 


GRAV. DIFF.: 

WEIGHT 

l 

HIG TO 

I. I. 

28.80 

1.0 

2 

HIG TO 

I. I. 

25.79 

1.0 

3 

BIG TO 

I. I. 

25.83 

1.0 

4 

BIG TO 

B.M. 

8.87 

1.0 

3 

HIG TO 

W/M 

S . 22 

1.0 

6 

HIG TO 

W/H 

5.23 

1.0 

7 

HIG TO 

*47 

-7.89 

1.0 

8 

BIG TO 

*171 

3.48 

1.0 

9 

HIC TO 

*324 

20.04 

1.0 

to 

HIG TO 

*323 

27.66 

1.0 

11 

*47 TO 

*171 

11.36 

1.0 

12 

*47 TO 

*171 

11.40 

1.0 

13 

*47 TO 

W/M 

13. 12 

1.0 

14 

*47 TO 

W/M 

13. 12 

1.0 

15 

*47 TO 

B.M. 

13.74 

1.0 

16 

*171 TO 

W/M 

1.73 

1.0 

17 

*171 TO 

*325 

24. 14 

1.0 

18 

*171 TO 

*324 

16.54 

1.0 

19 

*171 TO 

*324 

16.31 

1.0 

20 

*325 TO 

B. II. 

-21.74 

1.0 

21 

*325 TO 

I. I. 

-1.05 

1.0 

22 

*323 TO 

HICK 

-0.43 

1.0 

23 

*325 TO 

*324 

-7.60 

1.0 

24 

B.M. TO 

W/M 

-0.63 

1.0 

25 

B.M. TO 

I. I . 

19.92 

1.0 

26 

I. I. TO 

W/M 

-20.57 

1.0 

27 

1. 1. TO 

HICK 

1.40 

1.0 

28 

I . I . TO 

HICK 

1.37 

1.0 

29 

I. I. TO 

HICK 

1.37 

1.0 

30 

B.M. TO 

HICK 

21.32 

1.0 

31 

*171 TO 

HICK 

23.73 

1.0 

32 

*171 TO 

C.M. 

2.38 

1.0 

33 

B.M. TO 

*324 

14. 14 

1.0 

34 

*324 TO 

I. I. 

5.75 

1.0 

35 

HICK TO 

W/M 

-21.97 

1.0 

36 

*171 TO 

I. I. 

22.27 

1.0 

37 

*324 TO 

HICK 

7. 12 

1.0 

38 

HICK TO 

*47 

-35 . 06 

1.0 

39 

*47 TO 

I. I. 

33.69 

1.0 

40 

HICK TO 

HIG 

-27.21 

1.0 

CO! 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
17 

7DITI0N RESIDUALS BEFORE ADJUSTMENT: 
-0.03 
-0.03 
0.03 
0.02 
0.01 
0.00 
0.04 
0.00 
0.03 
0.03 
-0.03 
0.00 
0.00 
-0.03 
0.01 
-0.03 
0.00 


42 


18: -0.06 
19t-0.02 
20 : 0.00 
21: 0.08 
221-0.0S 

22i 0.00 
241-0.02 
23: 0.04 
26: 0.02 
27: 0.02 
281-0.01 
29: 0.11 
80: 0.00 
31: 0.00 


32: 0.02 
CORRELATES: 


l: 

0.016 

2: 

0.002 

3: 

-0.017 

4: 

-0.012 

3: 

0.038 

6: 

-0.008 

7: 

-0.039 

8: 

0.030 

9: 

0.017 

10: 

-0.050 

11: 

0.020 

12: 

0.005 

13: 

0.005 

14: 

-0.006 

15: 

-0.016 

16: 

0.014 

17: 

-0.011 

18: 

0.019 

19: 

-0.021 

20: 

0.014 

21: 

-0.011 

22: 

0.034 

23: 

-0.006 

24: 

-0.022 

25: 

-0.015 

26: 

-0.002 

27: 

-0.006 

28: 

0.024 

29: 

-0.055 

30: 

0.004 

31: 

0.041 

32: 

-0.013 

CORRECTIONS 

TO LIN1 


l: 

0.002 

2: 

0.012 

3: 

-0.028 

4: 

0.008 

5: 

0.016 

6: 

-0.014 

7 : 

0.011 

8: 

0.034 

9: 

-0.006 

10: 

-0.022 

11: 

0.019 

12: 

-0.021 
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13: -0.005 

14: -0.005 

15: 0.017 

10 

0.006 

17 

-0.002 

18 

-0.006 

10 

0.024 

20 

-0.020 

21 

0.014 

22 

-0.013 

23 

-0.004 

24 

-0.012 

25 

0.004 

26 

0.004 

27 

-0.007 

28 

0.023 

20 

0.023 

30 

-0.003 

31 

-0.085 

32 

-0.002 

33 

0.016 

34 

0.018 

35 

0.011 

36 

0.032 

37 

0.041 

38 

-0.014 

39 

-0.009 

40 

0.015 


ORIGINAL PAGE IS 
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sun w= 0.014 
check on sun w* 0.014 

HU* 0.021 


CORRECTED GRAVITY DIFFERENCES: 


* 

LINE: 



GRAV. DIFF.: 

WEIGHT 

1 

HIG 

TO 

I. I. 

25.80 

1.0 

2 

HTC 

TO 

I. I. 

25.80 

1.0 

3 

HIG 

TO 

I. I. 

25.80 

1.0 

4 

HIG 

TO 

B.M. 

5.88 

1.0 

5 

HIG 

TO 

W/M 

5.24 

1.0 

6 

HIC 

TO 

W/M 

5.24 

1.0 

7 

HIG 

TO 

*47 

-7.88 

1.0 

8 

HIG 

TO 

*171 

3.51 

1.0 

9 

HIG 

TO 

*324 

20.03 

1.0 

to 

H'C 

TO 

*325 

27.64 

1.0 

11 

*47 

TO 

*171 

11.38 

1.0 

12 

*47 

TO 

*171 

11.38 

1.0 

13 

*47 

TO 

W/'M 

13. 1 1 

1.0 

14 

*47 

TO 

W/'M 

13. 11 

1.0 

15 

*47 

TO 

B.M. 

13.76 

1.0 

16 

*171 

TO 

W/’M 

1.74 

1.0 

17 

*171 

TO 

*325 

24. 14 

1.0 

18 

*171 

TO 

*324 

16.53 

1.0 

19 

*171 

TO 

*324 

16.53 

1.0 

20 

*325 

TO 

B.M. 

-21.76 

1.0 

21 

*325 

TO 

I. I . 

-1.84 

1.0 

22 

*325 

TO 

HICK 

-0.44 

1.0 

23 

*325 

TO 

*324 

-7.60 

1.0 

24 

B.M. 

TO 

W/M 

-0.64 

1.0 

25 

B.M. 

TO 

I. I. 

19.92 

1.0 

26 

I. I. 

TO 

W/M 

-20.57 

1 .0 

27 

I. I. 

TO 

HICK 

1.39 

1.0 

28 

I. I. 

TO 

HICK 

1.39 

1 .0 

29 

I. I. 

TO 

HICK 

1.39 

1.0 


3# 

B.M. 

TO 

HICK 

21.32 

1.0 

31 

*171 

TO 

HICK 

23.79 

1.0 

32 

*171 

TO 

B.M. 

2.38 

1.0 

33 

B.M. 

TO 

*324 

14. 16 

1.0 

34 

*324 

TO 

I. I. 

3.77 

1.0 

33 

HICK 

TO 

W/M 

-21.96 

1.0 

36 

*171 

TO 

I. I. 

22.30 

1.0 

3T 

*324 

TO 

HICK 

7. 16 

1.0 

38 

HICK 

TO 

-47 

-33.07 

1.0 

34 

*47 

TO 

I. I. 

33.68 

1.0 

40 

HICK 

TO 

HIC 

-27.20 

1.0 


CONDITION RESIDUALS AFTER ADJUSTMENT: 


1 

: 0 

.00 

2 

: 0 

.00 

3 

: 0 

.00 

4 

0 

.00 

3 

0 

.00 

6 

0 

.00 

7 

0 

.00 

8 

0 

.00 

9 

0 

00 

10 

0 

.00 

11 

0 

00 

12 

0 . 

00 

13 

0 . 

00 

14 

0 . 

00 

13 

0 . 

00 

16 

0 . 

00 

17 

0 . 

00 

18 

0 . 

00 

19 

0 . 

00 

20 

0 . 

00 

21 

0 . 

00 

22 

0 . 

01 

23 

0 . 

00 

24 

0 . 

00 

23 

0 . 

00 

26 

0 . 

00 

27: 

0 . 

00 

28: 

0 . 

00 

29: 

0 . 

00 

30: 

0 . 

00 

31 : 

0 . 

00 

32: 

0 . 

00 
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ABSTRACT 


This report describes the planning, observations and 
adjustment of high-precision gravity survey networks established 
on the islands of Maui and Oahu as part of the geodetic- 
geophysical program in support of lunar laser ranging at 
Haleakala, Maui, Hawaii. 

The gravity survey networks include 43 independently 
measured gravity differences along the gravity calibration 
line from Kahului Airport to the summit of Mt. Haleakala, 
together with some key points close to tidal gauges on Maui, 
and 40 gravity differences within metropolitan Honolulu on 
Oahu . 


The results of our 1976-1978 survey are compared with 
surveys made in 1961 by the 1381st Geodetic Survey Squadron, 

Air Photographic and Charting Service, United States Air 
Force and with those made in 1964-1965 by the personnel of 
Hawaii Institute of Geophysics on the islands of Maui and Oahu. 

All final gravity values are given in the system of the 
International Gravity Standardization Net 1971 (IGSN 71); 
values are obtained by subtracting 14.57 mga 1 from the Potsdam 
value at the gravity base station at the Hickara Air Force Base, 
Honolulu. For the annual comparison of inter-island gravity 
surveys the secondary gravity base station at the Kahului 
Airport, Maui has been selected as the reference station. The 
tie between the two gravity networks was made from Kahului 
Airport to the inter-island terminal at Honolulu International 
Airport . 
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I. INTRODUCTION 


( To verify and measure any tectonic plate motion between 

the Hawaiian Islands and the western United States mainland 
| and to collect data for polar motion studies, the United States 

lunar laser ranging program maintains two permanent lunar laser 
| ranging installations. The McDonald Observatory, on the North 

| American tectonic plate, is located at the summit of Mt. Locke 

I in the Davis Mountains of western Texas. It is utilised by 

the Marine Science Institute, University of Texas. LURE 
* Observatory, on the Pacific plate, is located at the summit 

| of Mt . Haleakala on the Island of Maui. It is utilized by 

5 the Institute for Astronomy, University of Hawaii. 

i To obtain maximum information about plate motion from 

lunar laser ranging In a volcanically active region such as 
the Hawaiian Islands, any local or regional motions should be 
determined. Through a joint effort of the University of Hawaii's 
f Institute of Geophysics and Institute for Astronomy, such 

motions are monitored through repeated geodetic laser surveys 
t between the LURE Observatory and selected points on Maui and 

neighboring islands, repeated first order levelings of key 
points on Maui, and repeated gravimetric surveys on Maui and 
f Oahu. Also used are ocean tide gauges, seismic monitoring of 

crustal activity, and tilt meter monitoring of changes in the 
\ local vertical (Carter et al., 1977). 

This report describes the results of gravimetric measure- 
f ments made in 1976-1978, compares the values with older measure- 

ments, and analyzes the reliability of the results. If the 
above procedure is done yearly, any significant changes in 
absolute gravity values could be the signal for further inves- 
tigations. Utilization of periodic first order leveling reveals 
whether changes in the absolute gravity value of a station are 
caused by a change in elevation or a change in the distribution 
* of masses beneath the station. In the case of no measurable 

; change in elevation it would be assumed that the distribution 

of masses had changed, and other geophysical investigations 
would be called for. Since yearly first order leveling would 
’ be costly, the gravity survey can also be used as a rough 

leveling survey to interpolate the periodic first order 
l leveling results. 


I 
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II. INSTRUMENTATION 

Gravity measurements were made with either o£ two gravity 
meters. The first was a La Coste and Romberg Model G Geodetic 
Gravity Meter No. 1; later measurements were made with a La 
Coste and Romberg Model G Geodetic Gravity Meter No. 144. The 
two gravity meters are similar in construction. Both have a 
reading precision of ± 0.01 milligal. The meters have a weight 
on the end of a horizontal beam supported by a zero-length 
spring. Their measuring system uses a measuring screw and 
lever system. Both meters have an internal heating unit to 
keep the temperature of the measuring system constant at 52.2°C 
for the G-l and 51.85 e C for the G-144. 

Operating procedure is the same for both meters. After 
the instrument is leveled, the measuring screw is turned to 
bring the crosshair viewed in the telescope eyepiece to the 
same spot on the scale for each reading. This means that for 
each measurement the weight is nulled at the same point. The 
reading taken is then a measure of the change of length of 
the spring between gravity points, and, therefore, of the 
change in gravity acting upon the weight. The correct calibra- 
tion factors must be applied to the dial reading obtained at 
each point, before earth tide corrections and instrument drift 
correction are applied. After these corrections have been 
made differential gravity values in milligals can be found 
between each measured point. 

The theoretical earth tide corrections were computed using 
a computer program assembled by H. C. Marsh of the Hawaii 
Institute of Geophysics in 1973. The program utilized the 
Longman equations written for computers. The vertical component 
of the lunar tidal force is given as (Woollard et al,, 1973): 


GM 


uMr 

7 ~ 


(3 cos 0 - 1) + 




Mr 


(5 c os 9-3 


cos 


9) (1) 


where (j, is Newton's gravitational constant, M is the mass of 
the moon, r is the distance from the observation site to the 
center of the earth, d is the distance from the center of the 
earth to the center of the moon, and 9 is the zenith angle of 
the moon. The vertical component of the solar tidal force is 
given as : 


GS * (cos 2 (|> - 1 ) 

D 


( 2 ) 


where S is the mass of the sun, D is 
center of the earth to the center of 
angle of the sun, and (j, and r are as 
component is, therefore, given as: 


the distance from the 
the sun, <p is the zenith 
above. The total vertical 
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GO - GM + GS (3) 

To Investigate the drift pattern of the La Coste and 
Romberg Model G Geodetic Gravity Meter No, 1 two sets of tests 
were made. During the first test period the Instrument was 
placed on the floor of Room 109A in the Hawaii Institute of 
Geophysics building and was carefully leveled. Before the 
start of the test the instrument had been kept continuously 
on heat for about one week. The test period lasted 174 hours 
and five to seven readings were made dally after careful 
leveling. The clamp of the horizontal beam was kept open 
continuously during the entire test period, the room tempera- 
ture was kept constant at 23°C during the day and at about 
25°C at night, and the temperature of the measuring system 
was a constant 52.2°C. In Table la, dates, time and number 
of observations are given, together with the daily mean and 
the maximum daily deviation from the mean. 

Twenty-one days after the end of the first set of tests 
the instrument was placed on the same location as before and 
identical testing procedure was repeated lasting this time 
295 hours. The results are given in Table lb. Between the 
two sets of tests the gravity meter was transported to Maul 
via air, used for seven days and transported around the island 
by car on very rough roads, and finally returned to Oahu via 
air. In spite of the rather small combined travel and non 
travel drifts shown in Tables la and lb we decided to close 
every survey loop to its initial point within one day. As an 
additional precaution, before opening the clamp of the horizon- 
tal beam after the gravity meter was properly leveled at gravity 
points, we turned the measuring screw to a precomputed or 
estimated correct position to prevent the crosshair bouncing 
within its mechanically limited tolerance. 
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TABLE 

la. Drift 

Pattern of G-l 

Gravity Meter 


Date of 
Observations 

Tine 

Mo. of Dally Mean, 

Observations mgal 

Maximum 

Daily 

Difference , 
meal 

7 .22-76 

0900-1600 

7 

2132.016 

+ 0.014 

7.23-76 

0900-1300 

5 

2.032 

- 0.012 

7.26-76 

1000-1630 

5 

2.094 

- 0.024 

7.27-76 

1200-1600 

5 

2.056 

- 0.005 

7.28-76 

1100-1545 

5 

2.064 

+ 0.016 

7.29-76 

1000-1500 

5 

2.068 

+ 0.012 


TABLE 

lb . Drift 

Pattern of G-l 

Gravity Meter 


Date of 
Observations 

Time 

No . of 

Observa t i ons 

Daily Mean, 
mga 1 

Maximum 

Daily 

Difference, 
mga 1 

8.19-76 

0900-1700 

9 

2131.890 

+ 0.010 

8.23-76 

1000-1500 

5 

1 .992 

+ 0,008 

8.24-76 

1000-1600 

6 

1.953 

- 0.013 

8.25-76 

0800-1600 

7 

1 . 944 

+ 0.016 

8.26-76 

1100- 1500 

4 

1 . 940 

- 0.010 

8.30-76 

0900-1600 

7 

1.973 

+ 0.017 

8.31-76 

0900-1600 

5 

1.980 

- 0.010 
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III. GRAVITY NETWORKS 

The network on Maul consists of nine gravity points selected 
so that It Includes five stations along the gravity calibration 
line from Kahulul Airport to the summit of Mt. Haleakala. These 
are: Station 1 (Kahulul Airport secondary base station)* station 

3* station 5* station IS (Haleakala National Park Headquarters), 
station 21 (Mees Observatory). The calibration line was estab- 
lished on Maul in 1964-1965 by the personnel of the Hawaii 
Institute of Geophysics to be used in testing the calibration 
constant of gravity meters. In addition, observations were 
made at new points in Hana Bay, La Perouse Bay, and Pier 2, all 
of them located next to tidal gauges. An extra gravity point 
was added to Hana Airport to be compared with surveys made at the 
same site in 1961 by the 1381st Geodetic Survey Squadron, Air 
Photographic and Charting Service, Unired States Air Force. 
Location of these gravity points can be seen in fig. 1, and 
further descriptions follow. 

1. Gravity station 1 at Kahulul Airport is in the terminal 
building along the wall near the airport manager's office. The 
station is marked by a small brass disk on the floor. 

2. Pier 2 station is at the western end of Pier 2 next to 
the USCGS bench mark #8 near the Kahulul tidal gauge. It is 
marked by a Phillip's head screw in a lead insert drll.ed into 
the concrete. 

3. Gravity station 3, located on the north side of Highway 
37 approximately 0.3 miles southeast of milepost 3, is marked 

by a 12-inch-square stepping-stone set in concrete at ground 
level . 


4. Gravity station 5 at the junction of Highways 37 and 
377 is in the middle of a triangular section of ground outlined 
by the junction of the two highways, and is marked by a 12-inch- 
square concrete block set at ground level. 

5. Station 15 at the Haleakala National Park Headquarters 
is on a stone curb approximately 20 m southwest of the entrance 
to the Park Headquarters building. Station is marked by an X 
chiseled into the stone curb. 

6. Station 21 in Mees Solar Observatory is located in the 
generator room of the solar observatory. Station is marked by 
a brass disk on the floor. 

7. Hana Bay station is on a concrete wall that leads toward 
the pier at Hana Bay. It is marked by a Phillip's head screw 

in a lead insert drilled into the concrete. The station coincides 
with the Hana Bay laser terminal, and is near the tidal gauge. 
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8. Has* Airport station it locatsd at tha southeast corner 
of an open lanai on the west side of the terminal building. 

9. La Perouse Bay station, at the end of Kihei Road by 
La Perouse Bay, is the southernmost of the three USCGS bench- 
narks located near the shoreline, and is near the tidal gauge. 

On Oahu, a nine-station network for gravity aeasureaenta 
was set up. All of the points were in the Honolulu area, fron 
Hickaa AFB to St. Louis High School. 

1. The Hickaa AFB station is inside the Military Air 
Coanand terainal building on the floor next to the doors leading 
from the U. S. Customs inspection arec. 

2. The inter-island terminal station Is located at the 
Honolulu International Airport beside a stone pillar at the 
inter-island terminal building. The point is narked by a brass 
disk on the sidewalk. 

3. Station UH 325 is on the southernmost part of Sand 
Island on a concrete pier at the shoreline. 

4. Station UH 324 is south of the road leading to the 
State park on Sand Island. The point is on a concrete slab 
near a wire fence approximately 20 m from the park entrance. 

5. The Bishop Museum station, in Room 2 of the Bishop 
Museum annex, is marked by a brass disk on the floor. 

6. The Office of Weights and Measures station is in the 
Hawaii State Department of Agriculture Weights and Measures 
building at 1428 S. King Street. The point la on the floor 
beside a large generator in the east room of the building. 

7. The Hawaii Institute of Geophysics station, outside 
HIG Room 108 near the wall, is marked by a brass disk on the 
sidewalk. 


8. Station UH 47 is located on the north side of Lowrey 
Street directly across from a church. The point is on a 12-inch- 
square concrete block next tc the curb at ground level. 

9. Station UH 171, on the west side of the road leading 
into St. Louis High School, is on a concrete block next to the 
curb at ground level approximately 50 m from Waialae Avenue. 
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IV. GRAVITY NETWORK ADJUSTMENTS 


The procedure utilised ie the ease for each gravity net- 
work. After measurements between the gravity stations have 
been made and corrected for earth tides and Instrument drl.t, 
a value for the gravity difference between stations can be 
calculated. These differences can be used to form loop condi- 
tions for the network, which can then be put into a least 
squares conditional adjustment process. After the adjustment 
process, a discrepancy-free network is obtained along with the 
most probable values of the gravity differences between stations 
If one of the station is assigned as a base along with an 
associated observed gravity value, the observed gravity value 
of each point can then be calculated. 

After all measured lines have been tabulated, the ccndition 
equatiors for the adjustment can be found. Table 2 gives each 
measured line, the date the measurement was made, the gravity 
difference between the terminal points In milllgals, and the 
Instrument used for the Maui gravity • t-aork (Schenck, 1978). 
Table 3 gives the same data for the (h'iu gravity network 
(Schenck, 1978). Figs. 2 and 3 show tue Maul and Oahu gravity 
networks, respectively, along with the number of measurements 
of each line. 

When developing the condition equations for the adjustment, 
we found it easiest to work with only triangular loops. After 
all measured lines have been used as one side of a loop, the 
number of conditions can be found easily by looking at the 
number of extra times each of the lines of the loop has beer 
measured. It is most important to remember that only lines 
that have actually been measured can be used when setting up 
the condition equations. For example, the loop:statlon 3 - 
station 21 - Hana Bay - station 3 was observed in the following 
way (Fig. 2). Station 3 was used as the base for the loop; the 
gravity meter was set up and readings were made at that station. 
Then the gravity meter was transported by car all the way to 
the summit of Mt. Haleakala to station 21 where it was set up 
and readings made. Since there are no direct roads between the 
top of Halecksla and Hana Bay, the gravity meter was transported 
back to station 3, bypassing it and proceeding along the Hana 
Highway to Hana Bay where readings were made. The gravity meter 
was then transported back to station 3 where readings were made 
and the loop was closed. In this manner the gravity difference 
between station 21 and Hana Bay became Independently observed. 

Another Important fact to remember is that each of the 
conditions must contain at least one measurement that does not 
occur in any of the other conditions. From the 43 measured 
lines of the Maul network, there were 35 condition equations that 
could be formed. The Oahu network contained 40 measured lines 
and 32 conditions. 
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TABLE 2, Maui Gravity Lines 


Date 

Line 

Gravity Diff, 

Instrument 

Line No 

8-76 

1-3 

+ 

27.45 

G-l 

1 

8-76 

1-3 

+ 

27.44 

G-l 

2 

9-76 

1-3 


27.45 

G-l 

3 

9-76 

1-3 

+ 

27.46 

G-l 

4 

5-77 

1-3 

+ 

27.45 

G-144 

5 

5-78 

1-3 

+ 

27.44 

G-l 

6 

5-78 

1-3 

+ 

27.44 

G-144 

7 

8-76 

3-5 

+ 

68.52 

G-l 

8 

8-76 

3-5 

+ 

68.53 

G-l 

9 

9-76 

3-5 

+ 

68.53 

G-l 

10 

9-76 

3-5 

+ 

68.54 

G-l 

11 

5-75 

3-5 

+ 

68.52 

G-l 

12 

5-78 

3-5 

+ 

68.55 

G-144 

13 

5-77 

5-15 

+ 

321.92 

G-144 

14 

5-78 

5-15 

+ 

321.89 

G-l 

15 

5-78 

5-15 

+ 

321.89 

G-144 

16 

5-77 

15-21 

+ 

240.66 

G-144 

17 

5-78 

15-21 

+ 

240.65 

G-l 

18 

5-78 

15-21 

+ 

240.68 

G-144 

19 

5-78 

3-H.B. 

- 

78.92 

G-144 

20 

8-76 

H.B.-H.A.P. 

+ 

9.98 

G-l 

21 

8-76 

3-H.A.P. 

_ 

69.01 

G-l 

22 

5-78 

3-H.A.P. 

- 

68.99 

G-144 

23 

8-76 

3-L.P. 

_ 

37.45 

G-l 

24 

8-76 

3-L.P. 

- 

37.46 

G-l 

25 

9-76 

1-L.P. 

_ 

9.98 

G-l 

26 

3-77 

1-L.P. 

- 

9.97 

G-l 

27 

5-78 

5-L.P. 

_ 

105.98 

G-l 

28 

5-78 

5-L.P. 

- 

105.98 

G-144 

29 

5-78 

15-L.P. 

_ 

427.87 

G-l 

30 

5-78 

15-L.P. 

- 

427.91 

G-144 

31 

5-78 

2-21 

+ 

663.70 

G-l 

32 

5-78 

3-15 

+ 

390.41 

G-l 

33 

5-78 

2-3 

+ 

32.64 

G-l 

34 


r 
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TABLE 2, Maul Gravity Lines (Cont.) 


Date 

Line 

Gravity Dlff. 

Instrument 

Line No 

5-78 

2-L.P. 


4.84 

G-144 

35 

5-78 

2-L.P. 

- 

4.83 

G-144 

36 

5-78 

21-L.P . 

- 

668.58 

G-144 

37 

5-78 

2-5 

+ 

101.15 

G-l 

38 

5-78 

1-2 

- 

5.16 

G-144 

39 

5-78 

3-21 

+ 

631.08 

G-144 

40 

5-78 

3-21 

+ 

631.10 

G-144 

41 

5-78 

21-H.B. 

— 

710.00 

G-144 

42 


5-78 


21-H.A.P 


700.09 


G-144 


43 
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TABLE 3. Oahu Gravity Lines 


Date 

Line 

Gravity Dlff. 

Instruaent 

Line 3 

11-77 

HIG-II 

+ 

25.80 

G-l 

1 

11-7/ 

HIG-II 

+ 

25.79 

G-l 

2 

12-77 

HIG-II 

+ 

25.83 

G-l 

3 

12-77 

HIG-BM 

+ 

5.87 

G-l 

4 

11-77 

HIG-W/M 

+ 

5.22 

G-l 

5 

12-77 

HIG-W/M 

+ 

5.25 

G-l 

6 

12-77 

HIG-M7 

- 

7.89 

G-l 

7 

12-77 

HIG-#171 

+ 

3.48 

G-l 

8 

12-77 

HIG-#324 

+ 

20.04 

G-l 

9 

12-77 

HIG-#325 

+ 

27.66 

G-l 

10 

12-77 

#47-#171 

+ 

11.36 

G-l 

11 

5-78 

#47— #171 

+ 

11.40 

G-144 

12 

12-29 

#47-W/M 

+ 

13.12 

G-l 

13 

5-78 

#47-W/M 

+ 

13.12 

G-144 

14 

2-78 

#47-BM 

+ 

13.74 

G-l 

15 

12-77 

#171-W/M 

+ 

1.73 

G-l 

16 

12-77 

#1 71—# 325 

+ 

24.14 

G-l 

17 

12-77 

#1 71— #324 

+ 

16.54 

G-l 

18 

12-77 

#171-# 324 

+ 

16.51 

G-l 

19 

12-77 

#325-BM 

_ 

21.74 

G-l 

20 

12-77 

#325-11 

- 

1.85 

G-l 

21 

12-77 

#325-Hick. 

- 

0.43 

G-l 

22 

12-77 

#325-#324 

- 

7.60 

G-l 

23 

12-77 

BM-W/M 

_ 

0.63 

G-l 

24 

12-77 

BM-II 

+ 

19.92 

G-l 

25 

11-77 

II-W/M 

_ 

20.57 

G-l 

26 

12-77 

II-Hick. 

+ 

1.40 

G-l 

27 

5-78 

II-Hick. 

+ 

1.37 

G-144 

28 

5-78 

II-Hick. 

+ 

1.37 

G-144 

29 

2-78 

BM-Hick. 

+ 

21.32 

G-l 

30 

2-78 

#171-Hick. 

+ 

23.75 

G-l 

31 

5-78 

#171-bM 

+ 

2.38 

G-144 

32 

5-78 

BM-#324 

+ 

14.14 

G-144 

33 

5-78 

#324-11 

+ 

5.75 

G-144 

34 

5-78 

Hick.-W/M 


21.97 

G-144 

35 




TABLE 3. 

Oahu Gravity Lines 

(Cont . ) 


Date 

Line 

Gravity Diff. 

Instrument 

Line No. 

5-78 

#171-11 

+ 22.27 

G-144 

36 

5-78 

#324-Hick. 

+ 7.12 

G-144 

37 

5-78 

Hick.-#47 

- 35.06 

G-144 

38 

5-78 

#47-11 

+ 33.69 

G-144 

39 


5-78 


Hick.-HIG 


- 27.21 


G-144 


40 





■Ml 


Maui gravity network (from Schenck, 1978). 
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The adjustment procedure 
tial leveling network, except 
differences are equal to one. 
obtained as 


is similar to that of a differen 
that the weights of all gravity 
The number of conditions, r, is 


r ■ n - u (4) 

where n is the total number of observed gravity differences 
and u is the number of gravity differences needed to determine 
observed gravity values to all stations with reference to the 
selected base station. From Figs. 2 and 3 it is easy to see 
that in both gravity networks, eight gravity differences were 
needed for that determination. This is a good check to the 
method described above in determining the number of conditions 
in complex networks: On Maui, r ■ 43 - 8 ■ 35 and on Oahu, 

r - 40 - 8 - 32. 

From the loop (disclosures the condition equations are 
formed : 


/ a v + w 
^ j j 1 


l Vj 


w. 


0 

0 


(5) 


u TjV^ + w r = 0, j - 1,2,3, ,...,n 

where a^,bj ... are known constants (in this case +1, - 1, or 
0), w^,W£ ••• are (disclosures, and v^,V 2 ... are corrections 
that when applied to the observed gravity differences will make 
the misclosures zero. 


In Eq. 5, n is larger than r and the so-called conditional 

n 2 

minimum adjustment is applied. In this case 2 v. = minimum 

J-l J 

must be satisfied simultaneously with the condition equations 
in Eq. 5. If these equations are multiplied by i ■ 

1,2,3, ...,r it is obvious that the following expression is valid 



„ rr 

2K 1 L a j v j 




Vj 



• • • 


- 2K 


l 


) 

L 


r v 

j j 



minimum 


( 6 ) 


